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Abstract 

We consider possible interpretations of the recently detected X-ray afterglow from the 
gamma-ray burst source GRB 970228. Cosmological and Galactic models of gamma-ray 
bursts predict different flux and spectral evolution of X-ray afterglows. We show that 
models based on adiabatic expansion of relativistic forward shocks require very efficient 
particle energization or post-burst re-acceleration during the expansion. Cooling neutron 
star models predict a very distinctive spectral and flux evolution that can be tested in 
current X-ray data. 
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1. Introduction 

The discovery by SAX of an X-ray afterglow from GRB 970228 (Costa et al. 1997b, c, 
hereafter C97b,c) provides crucial information for the understanding of gamma-ray bursts 
(GRBs). This detection was possible because of a relatively rapid response in pointing to a 
GRB error box within a few hours. The results of two 'fast' Target of Opportunity (TOO) 
SAX pointings of GRB error boxes are currently available. An observation within 16 hrs 
of the ~ 15arcmin 2 error box of GRB 970111 (Hurley et al. 1997a) did not detect X-ray 
emission above a flux of 10~ 13 erg cm -2 s -1 (2-10 keV band) (Costa et al. 1997b). On the 
contrary, a pointing within 8 hours at the ~ 6 arcmin 2 error box of GRB 970228 (Hurley et 
al. 1997b) resulted in the remarkable discovery of transient X-ray emission from a previously 
unknown source (C97c, Boiler et al. 1997). The average X-ray flux level (2-10 keV band) of 
the first TOO pointing at the GRB 970228 error box is ~ 3 x 10 -12 erg cm -2 s . A second 
observation carried out ~3 days later detected the same X-ray source at a substantially 
lower flux near 10~ 13 erg cm -2 s _1 (C97c). In this Letter we consider theoretical models 
that can reproduce the observed decay in X-ray flux by a factor of ~ 30 between the first 
(at t ~ 8 hrs, with t the elapsed time after the burst trigger) and second (at t ~ 3.6 days 
~ 86.6 hrs) SAX observations of GRB 970228. It is also interesting to note the absenceQof 
a relatively strong X-ray source (at flux level near 10~ 12 erg cm -2 s _1 ) ~ 16 hrs after the 
GRB 970111 trigger. The peak intensities of GRB 970111 and GRB 970228 as detected 
by the GRBM instrument in the 60-600 keV band are comparable within a factor of ~ 2 
(C97a,b). We infer that temporal evolutions of X-ray afterglows from GRBs differ among 
bursts, and suggest that GRB models with no universal afterglow evolution are most likely. 



1 Sources (a) and (b) detected by SAX in the WFC error box of GRB 970111 (Butler et 
al. 1997) are outside the refined error box of in 't Zand et al. (1997). 
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2. Cosmological models 

Coalescence of compact stars, failed supernovae, new-born spinars or special jet 
phenomena in distant galaxies can give rise to a rapidly expanding relativistic fireball 
dissipating its energy in internal and external shocks. To avoid pair opacity effects, 
these models require a large value of the relativistic bulk Lorentz factor T ~ 10 2 — 10 3 
characterizing a fast-moving radiation front. The observed photon energies are upscattered 
by T compared to the ones produced in the comoving frame. The observed spectral energy 
flux F v is T 5 times its comoving value, a factor of T 2 from the increasing solid angle from the 
apparent transverse dimension of the fireball (r± = ctT), and a factor of T 3 from relativistic 
beaming (Rees 1966). Adiabatic expansion applies to fireball emission occurring within a 
time ~ 10 3 — 10 4 times longer than the burst T 90 duration (T 90 >, 20 s for GRB 970111 
and GRB 970228, C97a,b). The initial burst may be produced by 'internal' or 'external' 
shocks for an initial value of T = T a ~ 10 2 — 10 3 ( e.g., Meszaros & Rees 1997, hereafter 
MR97). For an external shock, the deceleration radius is r dec ~ (10 16 cm) (E^i/noY^T^ 2 ^ 3 , 
and the observed deceleration timescale t^ ec ~ r dec /(cT 2 ), where E 5 i is the initial 
fireball kinetic energy, and n a the particle number density in the surrounding medium 
(we use the convention A x = A/(10 x )). We expect a temporal dependence of the type 
T(t) ~ r~ 3 / 2 oc t~ 3 / 8 and r\\ = ctT 2 oc t 1 / 4 , with ry the longitudinal radius. If I' u i(t) is the 
(time-dependent) specific intensity in the shock comoving frame [y 1 = z//T), the observer 
frame flux at a frequency v is then F v oc r\T(tf I' u (t) ~ t 2 T(tf I' u (t) (MR97). The 
temporal evolution of I' v (t) depends on the adopted model of emission. 

The prompt burst high-energy behavior as observed in the 1 keV - 1 GeV is consistent 
with synchrotron emission of impulsively accelerated particles (Tavani 1996a,b, hereafter 
T96a,b), possibly upscattered by inverse Compton (IC) effects (e.g., MR97). We adopt here 
the synchrotron model of prompt GRB emission from a fireball, and follow the evolution 
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of the radiative shock as time progresses. IC is not expected to play an important role 
at the large values of rn implied by the X-ray afterglow of GRB 970228. A shock model 
of emission requires special conditions in order to have detectable X-rays 8-80 hrs after a 
GRB. Failed supernova (Woosley 1993), accretion-powered emission from massive black 
holes (e.g., Woosley 1996) and spindown-powered emission from newly born spinars (e.g., 
Usov 1992) may also produce delayed X-ray emission. 

Let us assume that most of the emission of an expanding and progressively decelerating 
fireball is produced by the reverse shock, i.e., by the ejecta particles cooled by adiabatic 
expansion. The comoving intensity at r\\ > rd ec turns out to be I' u , oc n' e j B' Ar y , where 
n' ej is the ejecta radiating particle density, B' the magnetic field, and Ar[| = r\\/Y. For 
adiabatic expansion, n' ej ~ V' 1 ~ t~ 9/8 , B' ~ V'~ 2/3 ~ t~ 3/4 , and Ar[| ~ t 5/8 , where we 
assumed a frozen- in magnetic field with V = 4nr'^ Arf| the comoving volume (MR97). 
We therefore obtain I' v , ~ r 5 / 4 , and F v ~ t^ s . The observed X -ray afterglow can be 
residual synchrotron emission of non-thermal particles accelerated by the prompt shock 
and adiabatically cooled. A non-thermal particle distribution function is assumed to be 
formed at the prompt GRB shock, as indicated by broad-band GRB spectra (T96a,b). 
This function is described by a relativistic Maxwellian low-energy component below the 
critical energy 7*, and by a power-law component up to a maximum energy 7 ma:r above 7*. 
The detectability of an X-ray afterglow from an expanding fireball of this kind crucially 
depends on the synchrotron critical frequency oc Y 7^ B' evolved at the radiating site. 
Since 7* ~ 1/r ~ t~ 1//4 , the critical frequency has a relatively fast decay, z/*(t) ~ t~ 13 / 8 . The 
observed X-ray flux f x (t) is the integral of F v (t) over a fixed energy band (say, 2-10 keV), 
and it turns out to be f x (t) = (, x (t) v*{t) F Ut (t), where £ x (t) gives the fraction of total 
luminosity in the required energy range. We obtain f x (t) = £ x (t)t~ 21 / 8 , i.e., a decay 
behavior substantially steeper than observed for GRB 970228 (but consistent with the lack 
of afterglow from GRB 970111). 
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Shock emission from an adiabatically cooled fireball is predicted to be non-thermal. 
However, due to the strong temporal dependence of 14, delayed X-ray emission from 
impulsively accelerated particles is possible only for very efficient energization. We can 
constrain the post-shock particle energization for no re- acceleration during adiabatic 
expansion at ry > r dec . The observed photon energy of the afterglow is assumed to 
be E max ~ r(t) 7^.(0 B'(t) ~ (t/t*)~ 13/8 - The decay timescale t* 

can be deduced 

from the time evolution of the peak energy of the v F v spectrum during a typical GRB, 
~ (200 keV) (t/t*)~ 13 / 8 . A typical value is t* ~ 3 s. For no re-acceleration, we obtain 

ImaM/lM Z 0.22(t/**) 13/16 (!) 

i.e., 7 m ax(0)/7*(0) ~ 400 for t = 8 hrs. A substantial fraction of the shock luminosity 
is required to be emitted in X-rays even for values of v„ substantially smaller than 
h ~ 1 keV. If we adopt F v ~ z/ -1 above (marginally consistent with the OSSE 
observation of GRB 970228 at t — 0.5 hrs, Matz et al. 1997), we expect the fraction £ x to be 
constant during fireball evolution as long as Eq. 1 is satisfied. A steeper F v above would 
result in an X-ray flux steeper than f x (t) oc t~ 21 / 8 . In this model, synchrotron emission 
peaks in the optical range at t/t* ~ 1.8 x 10 3 . However, due to the steep decrease of the 
bolometric luminosity, the delayed unabsorbed optical emission from pure synchrotron is 
too faint to be detectable. 

An alternative model of fireball expansion relies on the emission by newly shocked 
material of the surrounding medium. In this case, the density of the radiating particles of the 
comoving frame is n' = V n . In the absence of Rayleigh- Taylor instabilities at the contact 
discontinuity, the forward shock fluid does not mix with that of the reverse shock. The 
magnetic field can be deduced from equipartition arguments, B' ~ X^ 2 (8 ix n Y 7* m p c 2 ) 1 / 2 , 
with \b the efficiency of turbulent generation of magnetic field energy, and m p the proton's 
mass. It is interesting that in this case, I' v , ~ t~ 5 / 8 ,F v ~ t°, and ~ t -3 / 2 (e.g., MR97). 
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We therefore deduce the decay behavior of the bolometric synchrotron emission in the 
2-10 keV band as f x (t) = £, x (t)t~ 3 / 2 . For a relatively constant £ x (t), this function is 
marginally consistent with the observed afterglow of GRB 970228. As for the previous 
model, we can obtain the conditions on the acceleration mechanism that make possible 
delayed X-ray emission. If the local 7* ~ T, and •y max ~ a'T, with a 1 £ 1, we can translate 
the condition on the observed photon energy as a requirement on initial values of 7* and 
7 max . For E c = 10 keV and t = 8 hrs, we obtain T 4 ~ 4.7 x 10 18 , i.e, r o ~ 4.6 x 10 4 . More 
generally, we obtain for no re-acceleration 

7 m a,(0) 2 7 *(0) 1 / 2 * 4.7 x 10 15 n o -^ r -f X ~W {t/8 hrs) 3 / 2 (2) 

We deduce that a forward shock may be consistent with the observed afterglow of 
GRB 970228, but only if the particle acceleration process of the forward shock material is 
very efficient and/or r o is very large. Fig. 1 show the results of the calculated flux in the 
X-ray and optical R bands for two different choices of the power-law index of the post-shock 
particle distribution function corresponding to F u ~ z/ -1 and F v ~ z/~ 3 / 2 , respectively. 
From Fig. 1 we deduce that in the absence of additional acceleration, the predicted X-ray 
flux does not agree with observations of GRB 970228 if t* < 100 s. The resulting optical 
transient emission in the absence of absorption effects is shown in Fig. 1 (dot-dashed 
curves). The initial ratio of optical to X-ray flux is ^ 10 -5 , implying an optical transient of 
R-magnitude tur ~ 18.4. The transient evolves as t~ 3 ^ 2 up to t/t* ~ 3 x 10 3 corresponding 
to the time when sweeps the optical band. For later times, the evolution of the optical 
flux becomes steeper. Detectable delayed radio emission may also be expected in this class 
of models (Paczyhski & Rhoads 1993). 

If particle re-acceleration occurs during the adiabatic expansion phase, the constraints 
on r ) m ax{t) and 7*(t) become less stringent. Turbulent mixing of the reverse and forward 
shock fronts may favor efficient particle energization. In this case, Eq. |2| may not apply and 
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a milder constraint for r y max (t) and 7*(£) can be required. The evolution of both the X-ray 
and optical fluxes can be less steep than calculated here. 

Models based on relativistic jets from failed supernovae (Woosley 1993) are also based 
on forward shock emission. Depending on the density of the surrouding medium, dissipation 
by synchrotron (or bremsstrahlung) emission can result in different values of the decay 
constant t* compared to models considered above. However, the qualitative features of 
the X-ray evolution should follow the trend given above. We can deduce an interesting 
constraint on the 'spinar' model for GRBs. A newly-born strongly magnetized neutron star 
can spin down very rapidly and make possible an initial fireball and subsequent non-thermal 
'pulsar-like' emission (e.g., Usov 1992). In order to explain the required energy near 
10 52 ergs _1 , this model assumes extreme values of the surface magnetic field ~ 10 15 G and 
initial spin period P 1CT 3 s. The spindown energy varies with time as ~ (t/t sd )~ 2 , with 
t s d the spindown timescale, t S( i ~ P/P ~ 100 s. A fraction of the total spindown luminosity 
can be emitted as high-energy radiation. At t — 8 hrs, the emitted luminosity is decreased 
by a factor of ~ 10~ 5 compared to the initial one. At t = 85 hrs, the decrease factor is 
~ 10~ 7 . Therefore, this model predicts a luminosity evolution in disagreement with the 
observations of GRB 970228. 

3. Galactic models 

At a distance d = rfioo 100 kpc, an observed X-ray flux F x = F_ i2 10~ 12 ergcm~ 2 s" 1 
corresponds to an isotropically radiated luminosity, L x ~ 10 36 .F_i2 rf 2 00 erg s _1 . This 
sub-Eddington luminosity for a solar mass compact star can be produced in different ways. 
A sudden internal explosive event from a neutron star might cause matter liftoff and 
substantial release of energy in the external crust. Neutron stars in an extended Galactic 
halo require explosion energies of order of 10 41 — 10 42 ergs. A super-Eddington luminosity 
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leads to matter liftoff from the surface, possibly followed by matter fallback onto the surface 
of the neutron star. Let us consider first the case of internal energy release and cooling with 
no contribution from subsequent accretion. The thermal response of a neutron star crust 
depends on the depth at which most of the residual burst energy is liberated. For deep 
energy deposition, the timescale for conduct ion- driven cooling might be of orders of years 
or more. On the other hand, energy deposition in shallow layers of the external crust may 
lead to cooling timescales of order of 10 — 10 4 s. This can be seen by considering the initial 
photon diffusion time of the surface layers of a neutron star, Tdiff ~ 3 k px 2 u t h/ (c« 7 ), 
where k is the radiative opacity, p the mass density, x the depth (in cm), u t h the thermal 
energy, c the speed of light, and w 7 the radiation energy density (e.g., Eichler & Cheng 
1989). For 56 Fe matter, we obtain T diff ~ (30 s) (rr/100) 2 (T 8 /0.3)~ 5 - 5 p 2 - 4 . Thus a post-burst 
neutron star surface radiating near the Eddington limit can have T 8 ~ 0.3 and Tdiff ~ 30 s. 
Solutions of the conductive heat flow equation for vanishing depths depend on how the 
heat conductivity and heat capacity vary as a function of depth. This dependence is poorly 
known and we consider here the general form of the solution T(t) = T Q (t/V c )~ 1+ ^, where ( 
depends on the functional dependence of heat conductivity and capacity on depth. 

We can then assume that the post-burst bolometric flux from a neutron star surface 
as detected by a distant observer depends on time as L(t') = L Q t'~ a . This implies 
an effective surface temperature dependence T(t') = T Q t'~ a / 4 , where t' = t/r c and 
T oc L l J A (1 — i? S ch/-R*)~ 1//4 , with R* and R Sch the neutron star and Schwarzschild radii, 
respectively. For an initial luminosity radiated by the whole star near the Eddington limit 
Le = L 38 lO^ergs" 1 , we deduce an initial temperature, T e ~2x 10 7 L^ 8 4 K. It is interesting 
that for a decay timescale r c ~ 100 s and exponent a = 3/2, the temperature of a cooling 
neutron star surface would be in an optimal range for detection by X-ray instruments 
~1 day after the burst. Fig. 2 shows the calculated energy flux in the 2-10 keV energy 
band in units of the initial bolometric flux emitted by a cooling neutron star decaying as 
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fx(t') = fot'~ 3 ^ 2 , where f„ = L /(A-Kd 2 ). We consider two models characterized by different 
initial temperatures, T Q = 2 x 10 7 K and T Q = 3 x 10 7 K, respectively. We notice that the 
X-ray fluxes detected approximately 8 and 85 hours following GRB 970228 qualitatively 
agree with expectations based on this simple model. An interpretation in terms of thermal 
afterglow implies ^ — 5/8, a value that can be used for a detailed modelling of the thermal 
response of a neutron star. 

The cooling neutron star model predicts a substantial softening of the spectrum. Fig. 2 
also shows the calculated evolution of the hardness H defined as the ratio of photon fluxes 
in the 3-6 keV and 2-3 keV bands. The hardness ratio is predicted to vary by a factor ~ 2 
within the first SAX TOO observation of GRB 970228. An even more drastic variation of 
the hardness dropping below 0.1 is predicted for later times, e.g., during the second SAX 
TOO observation. No detectable optical emission is expected in this case from pure neutron 
star cooling. 

Other ways of producing X-ray emission can be considered. If the bulk of GRB 
emission occurs outside but not too far from a neutron star, photon irradiation and particle 
precipitation can further heat its surface. The energy deposition from above can lead to 
a thermal relaxation qualitatively similar to that of a cooling neutron star surface heated 
from below. However, the detailed heat transport may be different than in the previous 
case. Furthermore, fallback material from the burst explosion can settle in a disk or 
spherical inflow onto the neutron star surface. The X-ray spectrum in this case is predicted 
to be a combination of blackbody and power-law components as observed from several 
low-luminosity accreting neutron stars (White, Nagase & Parmar 1995). X-ray pulsations 
might be detectable in case of emission from a strongly magnetized neutron star. The ratio 
of optical to X-ray luminosity is expected to be 10 -2 — 10~ 3 , i.e., similar to faint accreting 
compact sources (e.g., White et al. 1995). 



4. Discussion and conclusions 



We showed that the detection of an X-ray afterglow from a GRB strongly constrains 
theoretical models. Table 1 summarizes the main properties of the models considered here. 
Cosmological models based on adiabatic evolution of a synchrotron radiating shock front 
are the most constrained. The calculated emission in models of adiabatically expanding 
reverse shocks does not agree with observations of GRB 970228. Forward shock models 
might be compatible with observations only for very efficient particle acceleration, large 
decay timescale t* >, 100 s or re-acceleration during expansion. An optical transient of 
initial magnitude m R ~ 18 — 19 lasting a few hours after the GRB is expected in this 
model. If t* = 300 s, we deduce from Fig. 1 an X-ray flux in qualitative agreement 
with the observations of GRB 970228 and an optical magnitude at t/t* ~ 100 equal to 
m R ~ 22.7, i.e., slightly higher than the limiting magnitude (m R ~ 22) of the observations 
at t — 15.3 hr reported by Guarnieri et al. (1997). Consistency with the optical transient 
possibly associated with GRB 970228 (Groot et al. 1997) requires t* >, 300 s. No radio 
transient source above ~1 mJy has been reported in searches at 5 GHz of the refined error 
boxes of GRB 970111 and GRB 970228 (Frail et al. 1997, Galama et al. 1997). We also 
showed that Galactic models of post-burst emission from cooling neutron stars have definite 
predictions. Different flux and spectral evolution patterns of X-ray afterglows are expected 
as a function of initial surface temperature and cooling timescale. 

Our analysis is important for the interpretation of follow-up X-ray observations of 
GRB 970228. As shown in Fig. 2, in neutron star cooling models no detectable X-ray flux 
above 10~ 14 erg cm -2 s -1 from GRB 970228 is expected in the X-ray band for t > 100 hrs. 
Analogously, cosmological models based on forward shocks with no re-acceleration, predict 
an evolution of the X-ray flux f x ~ t~ 3 / 2 or steeper. Again, the predicted flux at t >, 100 hrs 
is below detectability with current instruments. A violation of these predictions would 
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be remarkable, since it would point to a steady-state high-energy emission of the GRB 
counterpart. 

We emphasize that X-ray afterglows from GRBs show different time evolutions. 
Within a given model, variations of physical properties (critical timescales, temperature, 
maximum value of particle energies in non-thermal models, etc.) can explain the difference 
between the afterglows from GRB 970111 and GRB 970228. A combined sequence of fast 
multiwavelength observations following GRBs can further constrain the emission parameters 
and their variations among different GRB sources. 

The author thanks E. Costa, F. Frontera, M. Ruderman and M. Feroci for discussions. 
Research supported in parts by the NASA grant NAG5-2729. 
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FIG. 1 - Calculated evolution of X-ray and optical (R-band) fluxes for an adiabatically expanding 
forward shock with no re-acceleration. The temporal axis scale is in units of t*, with the 
synchrotron critical frequency evolving as ~ (t/t*r 3/2 . (Solid curve:) X-ray (2-10 keV) flux (in 
units of the initial bolometric flux) for a spectral form F v ~ v~ l above i/*; (long-dashed curve:) 
the same for F v ~ v~^l 2 above v*. (Long- dashed- dotted curve:) R-band flux (in units of the initial 
bolometric flux) for F u ~ above i/*; (short- dashed- dotted curve:) the same for 
above v*. 
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FIG. 2 - Evolution of X-ray flux and photon hardness as a function of tjr c for the cooling neutron 
star model. (Solid curve:) X-ray flux (2-10 keV) in units of F ^ \ for an initial temperature 
T Q = 2 x 10 7 K; (dashed- dotted curve:) the same for T Q = 3 x 10 7 K. (Short-dashed curve:) photon 
hardness ratio for the channels 3-6 and 2-3 keV for T Q = 2 x 10 7 K; (long-dashed curve:) photon 
hardness ratio for T Q = 3 x 10 7 K. 
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TABLE 1: Models for X-ray afterglows from GRBs 
Model fx{t) Requirements 



Adiabatic fireball: reverse shock 


£.r 21/8 


Adiabatic fireball: forward shock 




Failed supernova 




Spindown-powered emission 


r 2 


NS cooling 




NS cooling with external irradiation 


r i+C 



~ COnst.; 7maa;(0)/7*( ) ~ 400 

& ~ const.; 7m a:r (0) 2 7*(0) 1/2 > 4.7 x 10 15 n^ 2 T~ 3 2 

Optically thin circumstellar medium 

Strongly magnetized, rapidly rotating compact star 

Shallow energy deposition, ( ~ 5/8 

Substantial reprocessing, (' ~ 5/8 



